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We report here the discovery of a significant source of systematic error in the rotation measure 
determinations of pulsars. Conventional analysis of high sensitivity polarimetric observations of 
PSR B2016+28 display variation of the rotation measure of ±15 rad m~^ (around the mean value 
of -34.6 rad m~^) across the pulse profile. Analysis of single pulse data shows that this variation 
is an artifact of the incoherent superposition of quasi-orthogonal polarisation modes along with 
the frequency dependence of relative strength and/or quasi-orthogonality of the modes. Quasi- 
orthogonal polarization is common among pulsars, and therefore this effect needs to be taken 
into account in the interpretation of pulsar rotation measures. 

Subject headings: pulsars: polarisation - radiation mechanism: non-thermal - ISM: Magnetic field 
1. Introduction sity. e and rUe are the charge and the mass an 



ABSTRACT 



Pulsars play a major role in our understand- 
ing of interstellar magnetic field structure owing to 
measurements of Faraday rotation of the plane of 
linear polarization as a function of radio frequency 
introduced by the component of the field along the 
sight line in the warm interstellar medium 
The standard definition of rotation measure (RM) 
in c.g.s. units is given by 



electron, and c is the velocity of light in vacuum. 
The amount of rotation experienced by the intrin- 
sic linear polarisation position angle of the source 
(ipo) at a given wavelength (A) is expressed as 



(2) 



where L is the distance to the pulsar from the 
observer, dl is the length element along the line 
of sight, and rig is the free electron number den- 
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where ^p is the apparent linear polarisation posi- 
tion angle (PA) as seen by the observer. Several 
investigators in the past have measured RM values 
of pulsars (Manchester 1972; 1974; Hamilton et al. 
1981; Hamilton & Lyne 1987; Rand & Lyne 1994; 
Qiao et al. 1995; Weisberg et al. 2003). Time de- 
pendence in the measured RM values^ have also 



It is also important to compensate for the variable iono- 



been noted against some pulsars like the ones in 
the Vela and Crab supernova remnants (Hamilton 
et al. 1977; Rankin et al. 1988). A positive RM 
means the direction of magnetic field is towards 
the observer, and a negative RM means otherwise. 
Using the measured values, detailed modelling of 
magnetic field structure has been carried out by 
several authors in the past (Thompson & Nelson 
1980; Lyne et al. 1989; Indrani & Deshpande 1998; 
Han et al. 1999; 2002; Mitra et al. 2003). All of 
these galactic models based on pulsar RMs have 
a central assumption that the RM is completely 
determined by the interstellar medium, and that 
the magnetosphere of the pulsar, with all of its 
complexities, does not contribute significantly. 

We show in this paper that a conventional RM 
analysis based on average pulse profiles leads to 
large variations of the RM across the pulse profile 
of PSR B2016+28. The conventional analysis re- 
quires some form of averaging to arrive at a single 
RM. Additional work that we are doing has shown 
that this effect occurs in other pulsars as well. In 
pulsars with significant RM variations across the 
pulse profile, there is a source of error that may not 
have been considered in past modeling of galactic 
magnetic fields. 

Significant magneto-ionic propagation effects 
are not expected in pulsar magnetospheres owing 
to the ultra-rclativistic nature of the plasma. Any 
large Faraday rotation within the emission region 
of the magnetosphere would lead to severe dcpo- 
larisation across our band. The fact that pulsar 
radiation is highly polarised therefore shows that 
there is no significant Faraday rotation within the 
emission region. We arc led then to look more 
closely at the data to determine the nature of the 
RM variations across the pulse. 

We show in this work an improved method for 
RM determination when sufficient signal-to-noise 
ratio (S/N) allows detection of the pulsed radiation 
in single pulses. In this case, the orthogonal modes 
of polarization (e.g., Backer & Rankin 1980) can 
be identified, and the RM can be determined 
from each mode independently. If the modes were 
strictly orthogonal AND there was no jitter AND 

spheric contribution to the measured RM values. As 
Manchester (1972) states, this contribution could be as high 
as some 0.1 to 6 rad m~^, depending on the time of the 
day and the declination of the pulsar. 



there was no frequency dependence of these prop- 
erties, then this procedure would not be necessary. 
However, the PAs (position angles) at each longi- 
tude do display jitter and a slight degree of non- 
orthogonality (Gil, Snakowski & Stinebring 1991; 
Gil et al. 1992). We show in this work that the 
B2016-I-28 emission has a slight frequency depen- 
dence in the relative strengths and (or) amounts of 
non-orthogonality. This is the source of the pulse- 
longitude dependence of the RM in the conven- 
tional analysis. One does not have to appeal to 
strong, and unexpected, magneto-ionic propaga- 
tion effects. 

Most pulsars for which any amount of single 
pulse study exists show signs of orthogonal modes 
in their emission. The problem of origin of these 
modes has been addressed by several investigators 
(e.g., Melrose 1979; Petrova & Lyubarskii 2000; 
Radhakrishnan & Deshpande 2001; Petrova 2001). 
As they describe, the orthogonal mode could arise 
from the partial conversion of an original mode 
(ordinary) to the other (extra-ordinary) due to 
propagation in the magnetosphere of the pulsar. 
Even the cause for the slight non-orthogonality 
has been addressed by investigators like Petrova 
(2001). 

In this paper, after briefly describing the de- 
tails of our observations in Sec. 2, we summarize 
the conventional approach to RM determination 
in Sec. 3. Here we describe the aforementioned 
RM variations as a function of pulse longitude for 
PSR B2016-h28. In Sec 4 & 5 we show how quasi- 
orthogonal modes in pulsars can cause severe arti- 
facts in RM estimations, and in particular how the 
apparent variations in the RM of PSR B2016+28 
can be produced. We conclude our analysis with 
a detailed discussion in Sec. 6. The symbols for 
variables used in this paper are defined in the ap- 
pendix. 

2. Observations 

The average piilse observations were performed 
at the Arecibo Observatory at 430 MHz in May 
and December 1992 with a band width of 5 MHz. 
Average pulse profiles were produced by integrat- 
ing the signal in each frequency channel for 120 
seconds. The observation setup has been de- 
scribed in detail in Weisberg et al. (2003). 

The single-pulse observations were carried out 



2 



-10 -5 5 

longitude (deg) 



10 



-10 -5 5 

longitude (deg) 



10 




-5 5 10 

ongitude (deg) 



3 d 





1 1 1 1 1 1 1 1 1 1 1 




\ \ 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 








-10 



-5 

ongitude (deg) 



10 



Fig. 1. — (a) - (d) : Results from four different average-pulse data sets. Top panel gives the average profile 

in Stokes I parameter (solid line), linearly polarised power ('dash' line) and circularly polarised power ('dot- 
dash' line). Middle: Linear polarisation position angle determined from the average profile. Bottom: 
RM(0) vs 0. See text for details. 



at Arecibo Observatory at 430-MHz centre fre- 
quency in a single session on October 1992. With 
a special program for gating the 40-MHz corre- 
lator, auto-corrclation and cross-correlation func- 
tions of the right-hand and left-hand polarisation- 
channel voltages were recorded with 32 correla- 
tion lags. The bandwidth of the observation was 
10 MHz. The correlation vectors were averaged 
to a time resolution of 506 /xsec. The correla- 
tion data were thrcc-lcvcl corrected after scaling, 
and Fourier transformed to produce Stokes param- 
eter spectra. A detailed calibration procedure was 
adopted to correct for instrumental effects, delays 
and interstellar dispersion. These observations as 
well as calibration procedures have been described 
in Rankin & Rathnasrce (1995). The resulting po- 
larised pulse sequence had 32 channels (although 
channel 1 lacked Stokes U and thus was unusable) 
with 3043 pulses (1600 sec), and was gated to in- 
clude a longitude range of some 45°. A constant 
RM value of -34.6 rad m~^ (earlier measurement 
in the literature) was also subtracted from the 
data. 

3. Conventional RM measurements 

As individual pulses from pulsars are in general 
faint, one typically generates an "average" pulse 
profile by folding the time series in each frequency 
channel and Stokes parameter synchronously with 
the Doppler-shifted apparent period. An impor- 
tant property of radio emission from pulsars is that 
the position angle of linear polarisation changes as 
a function of pulse longitude (Radhakrishnan & 
Cooke 1969). Owing to this property, RM mea- 
surements of pulsars have always been more com- 
plex than for any other polarised radio source in 
the sky. PA measurements using a given longitude 
bin of folded profiles and all frequency channels 
are used to fit for the rotation measure [RM(0)] on 
the basis of Eqn. 2. This is repeated for each lon- 
gitude bin, and the RM from all these bins are 
averaged to estimate the overall value, (RM), for 
the pulsar. Of course, the underlying assumption 
here is that the value measured in all the longi- 
tude bins are entirely introduced by only the in- 
terstellar medium, and hence there is no longitude- 
dependent Rotation Measure introduced by the 
pulsar itself. Errors in the RM((/)) can be simply 
propagated using the S/M of each individual esti- 
mate and an assumption of independence of the 



estimates. 

In our analysis, as the first step, we have at- 
tempted to compute RM(0) values as a function 
of pulse longitude. As mentioned above, this 
has been estimated from the average profiles con- 
structed in each of the 31 frequency channels in 
our sequence. Fig. 1 & 2 summarise the re- 
sults. Fig. 1 shows the results from four differ- 
ent scans from our, each of 120 sec integration. 
In the top panel, we give average pulse profile in 
Stokes I parameter (solid line), average linearly 
polarised power ('dash' line), and average circu- 
larly polarised power ('dot-dash' line). The mid- 
dle panel shows linear polarisation position angle, 
and the bottom panel shows RK{^). 

In the top panel of Fig. 2, the solid line gives 
the average pulse profile in total power and the 
'thin' dash line gives the average profile in lin- 
ear polarisation computed from average profiles 
in all Stokes parameters. This is to be compared 
with the 'thick' dash line, which gives linear power 
computed from single-pulses and averaged over all 
pulses. As we can see, the 'thin' dash line shows a 
smaller degree of linear polarisation owing to in- 
coherent superposition of the polarisation modes. 
Then the 'thin' dot-dash line indicates power in 
circular polarisation computed from average pro- 
files in Stokes-V parameter. Power less than zero 
just means that it is left-circularly polarised. 

The dots in the middle panel give the linear 
polarisation PA curve as defined by the average 
profile, and the gray scale gives the probability 
density function (PDF) of the position angle com- 
puted from all the single-pulses. While comput- 
ing this PDF, we have weighted the values with 
the square of their S/N ratio of the polarised flux 
as defined by \/Q'^ + U^/asys (where Q k, U are 
two of the Stokes parameters, and (Jsys is the root- 
mean-square value of the system noise flux^). Al- 
though the PA range outside ±90° is redundant, 
we have chosen a range of ±180° for clarity. 

In the bottom panel of Fig. 2, we give the mea- 
sured values of RM(^) as a function of pulse longi- 
tude. The dotted line in the bottom panel shows 
the earlier measurement of RM for this pulsar 

■^As the telescope gain of Arecibo Telescope is significantly 
high, variation of system noise with the strength of the 
pulsar was explicitly taken into account with the help of 
calibrated average pulse profile. 
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Fig. 2. — RM(^) vs 4>. Top : Average pulse profile in Stokes-I is given as solid line. 'Thin' and 'thick' 
dash lines are linear polarisation power computed from average pulses and from single-pulses and averaged 
over all pulses, respectively. 'Thin' dot-dash line gives circular polarisation power (Stokes V) computed from 
average pulses. A horizontal dotted line is drawn to indicate the zcro-lcvcl in the Y-axis. Middle : The 
dots show position angle estimated from the average profile. The gray scale gives the probability density 
function of position angles estimated from single-pulses. Bottom : RM(0) estimated from average pulse. 
The value of (RM) computed from this is -38.3(3) rad m~^. 



(Manchester 1972), -34.6 rad m" 2. As we men- 
tioned in Section 2, the data that we used for 
our analysis already had a constant RM value of 
-34.6 rad m~^ removed from it. But to be con- 
sistent, we have added this offset in the plot (dot- 
ted line). Effective RM from our measurement 
comes to -38.3(3) rad m~^, which is significantly 
different from the earlier measurement of Manch- 
ester, as well as of Weisberg et al. (2003), which 
is -27.3 ±2.1 rad m'^. 

There are two aspects of Fig. 1 & 2 that needs to 

be discussed here. The measured values of RM((/)) 
as a function of pulse longitude display significant 
systematic variations. If this effect is true, then it 
has a very fundamental significance, as one docs 
not expect the interstellar medium to distinguish 



between one pulse longitude to the other. There- 
fore, this RM difference must be due to intrinsic 
propagation effects in the pulsar magnetosphcrc. 
In fact, this would become the first ever direct evi- 
dence for propagation effects in the pulsar magne- 
tosphere. However, as we will show in the follow- 
ing section, this is due to artifacts introduced by 
superposition of the two more or less orthogonal 
emission modes seen in this pulsar. 

Secondly, in Fig. 1, the overall behaviour of 
RM(0) is similar in all plots. As described in Sec. 2, 
these data sets were obtained in May and De- 
cember 1992. The "anti-symmetric" behaviour of 
RM(0) (with respect to longitude value of zero in 
Fig. 1) seems to be stable over a time scale of 
several months. However, although the overall be- 
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Fig. 3. — Rotation Measure of the two polarisation modes [TIA4a{4') ^ ^^.A^b (</')] as a function of pulse 
longitude. For comparison, we have given average profiles in total power (solid), linear (dash) and circular 
(dotted) polarisation in the top panel, and position angle PDF in the middle panel as gray scale. The bottom 
panel gives the Rotation Measure of the two modes, TZMa{4') and 72.A^b (</>)• Effective values of RM of the 
two modes are -35.0(4) & -34.5(3) rad uT^. 



haviour of RM(0) is stable, it is quite clear that 
the exact value of RM(0) is not the same in all the 
panels of Fig. 1. We will return to this later in 
Sec. 6. 

RM values varying as a function of pulse longi- 
tude is not unique to this pulsar alone. It turns 
out that this has been observed in other pulsars 
like PSR B0329+54 (Mitra, Johnston & Kramer, 
private communication) . 

4. Intrinsic RM((/)) variations? 

If the two polarisation modes are strictly or- 
thogonal, then they have no effect on conventional 
RM measurements. When we compute the aver- 
age profile in each frequency channel, all that we 
are doing is "incoherently" superposing the radia- 
tion from the two orthogonal modes. This would 



mean that the net degree of polarisation is the 
difference of their degrees of polarisation, and the 
PA direction is the same as that of the dominant 
mode. Provided adequate resultant polarised sig- 
nal, we can still estimate the RM value. If the 
modes are not strictly orthogonal, then the net 
PA of the average is the result of a vector sum- 
mation. This too is of no consequence to RM de- 
termination as long as the relative strength of the 
two modes and their individual mean PA remain 
constant as a function of radio frequency. 

There have been only a few quantitative stud- 
ies of the frequency dependence of the orthogo- 
nal modes (Karastergiou 2002; 2003). In general 
the degree of polarization of this pulsar decreases 
significantly with increasing frequency (Gould & 
Lyne 1998). Also, as the observations of Gould & 



6 



Lyne show, the mode-dominance transition point 
(longitude of ~ 29° in Fig. 2) moves to earher lon- 
gitudes with increasing radio frequency. This in- 
dicates that the relative strengths of the modes 
do not remain constant as a function of frequency. 
In other words, the first half of the profile, which 
is predominantly polarised with mode-A (repre- 
sented by a dash line) has a steeper spectral index 
when compared to the second half that is predom- 
inantly polarised with mode-B (dot-dash line). 

In the case of PSR B2016-I-28, it is also clear 
that the modes are not exactly orthogonal. The 
best fit to the two gray-scale tracks show that the 
PA and the slope at 30° longitude are : -42°.8(2) 
and 5.04(1)°/° for mode-A (represented by the 
dash line) and 49°. 8(1) and 4.37(1)°/° for mode-B 
(dot-dash line). These two slope values are signif- 
icantly different from each other. This PA non- 
orthogonality is known in the literature (Gil, 
Snakowski & Stinebring 1991; Gil et al. 1992), and 
has recently been studied in detail by McKinnon 
(2003). In fact, it has been shown by McKinnon 
that this is fairly common among several pulsars 
(e.g. PSRs B0950-I-08, B1929-M0, B2020-I-28). 

5. RM determination in presence of quasi- 
orthogonal modes 

Although the fractional polarisation seen in the 
average pulse profiles are typically a few tens of 
percent, individual pulses can exhibit even higher 
degree of polarisation, some times nearly reaching 
hundred percent (Stinebring et al. 1984; Rankin 
et al. 1989; Ramachandran et al. 2002; Rankin & 
Ramachandran 2003). Given the presence of the 
quasi-orthogonal modes and possible depolarisa- 
tion due to their superposition, quite independent 
of any possible artifacts, it makes sense for us to 
determine the RM with individual pulses rather 
than average pulses. 

As described in Section 2, our single pulse data 
set of PSR B2016-I-28 consists of 3843 pulses, 
gated to represent a longitude range of about 45° 
in each pulse. As the first step, in each frequency 
channel, we produced a PA distribution for every 
longitude. In doing this, we weighted each sam- 
ple by the square of S/N of the polarised fiux. At 
a given longitude, due to the presence of the two 
modes, one expects two significant peaks in this 
distribution separated by almost 90°. 
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Fig. 4. — Average deviation angle (Y-axis) from 
non-orthogonality of the two polarisation modes 
as a function of pulse longitude (X-axis). See text 
for details. 

As the second step, by a weighted mean, we 
found the ccntroid of the two peaks [xk{4'^ ^) ^^id 
Xb(0, J^)], both being functions of pulse longitude 
and radio frequency. We thus determined the 
mean PA values of each mode. 

Once the PAs are determined, then it is straight 
forward to fit for rotation measure values of the 
two modes, 7?,A^a(<^) and TZM.^{(j)), as a fimc- 
tion of pulse longitude. These two functions de- 
termined for PSR B2016-I-28 arc given in Fig. 3. 
Comparing this to the bottom panel of Fig. 2, it is 
obvious that the deviations of RM values observed 
for the two modes is far less than what is deter- 
mined with the average profiles. The effective RM 
values of the two modes turn out to be —35.0 ±0.4 
and —34.5 ±0.3 rad m~^, respectively. These val- 
ues arc very close to the earlier measiircment of 
-34.6 ± 1.4 rad m'^ (Manchester 1972). 

As emphasised in Section 4, the presence of the 
two modes itself, even if they are not orthogonal, 
should not introduce any artifact in RM measure- 
ments, as long as there is no frequency depen- 
dence of the relative strength and the intrinsic 
PAsof the two modes. However, in our analysis 
of PSR B2016±28, the value of RM(0) is signifi- 
cantly different from the values of 7?,A^a('/>) and 
7^A^b(0), which clearly indicates that the inco- 
herent superposition of the modes while generat- 
ing average profiles (as a function of frequency) 
has introduced a frequency dependence of the re- 
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sultant PA , which manifests as extra RM. 

To investigate this subtle effect in detail, with 
the help of the above described procedure, we pro- 
duced Fig. 4, that shows the non-orthogonality of 
the two modes as a function of pulse longitude. 
We have plotted the quantity 

= \- ([Xb(</', v) - Xa(</>, ■ (3) 

against pulse longitude (0). As wo can see, the 
deviations are not uniform across the pulse profile, 
and the variations are statistically significant. 

Let XpJ^^v) and X-Q{t,v) be the "instanta- 
neously" randomly varying linearly polarised in- 
tensities of the two quasi-orthogonal modes. Let 
C,{y) be the angle between the two vectors in the 
Stokes Q U space. Assuming that the PA of the 
primary mode is zero, we can write the two Stokes 
parameters Q and U as (see also McKinnon 2003) 

Q{t,v) = Xf^{t,u)+X^{t,v) cosC{y)+XQ{t,v) 
U{t,u) = XB{t,iy) smC{u)+Xv{t,iy) (4) 

If ({v) is independent of XA{t, v), then the PA as 
a function of frequency with average values of Q 
and U is 



(Xab(z^)) = ^tan" 
1 

= :rtan 



-FAB(i^) +C0SC(I^) 



(5) 



where the angular brackets indicate time-averaging, 

and Fab{i^) - {XA{t,iy))/{XB{t,iy)). X(^{t,u) 
and X\]{t,v) are the system noise strengths in Q 
and U Bs, & function of time and frequency. As we 
can see, in principle, frequency dependence in 
or i^AB(i') can introduce a frequency dependent 
(XAB(i^)), which can corrupt our rotation measure 
determination. 

The dependence of (xab(z^)) on frequency is dif- 
ferent from that of the RM. Therefore, in principle, 
the RM measured at different frequency ranges 
should be different. To check this, we divided our 
frequency band into two parts. Indeed, the value 
of RM((/)) was different in the two halves. In the 
first half of the band, the measured value of RM 
was —32.5(9) rad m~^, whereas in the other it was 
-40.5(4) rad m"^. The overall value of (RM) was 
-38.3(6) rad m'^. 



The error introduced by this effect depends on 
the bandwidth, centre frequency, frequency depen- 
dences of Q{v) and Fab(i')- Fig- 5 succinctly sum- 
marises this effect. This theoretically generated 
plot corresponds to a centre frequency of 430 MHz 
and a band width of 10 MHz, which is divided into 
128 frequency channels. In the left panel, we give 
C(i/) in the X-axis and effective RM introduced 
by this quasi-orthogonal mode-mixing in the Y- 
axis. The five curves, solid, dash, dot-dash, dot- 
ted and dot- dot- dot- dash, correspond to variations 
of fractional strength of the two modes [Fab(z^)] 
from the lower to the upper end of the frequency 
band of 0.3-0.9, 0.45-0.9, 0.6-0.9, 0.75-0.9 and a 
constant 0.9, respectively. As is seen, when the 
relative strength remains constant, then the RM 
introduced by this effect is zero. However, when 
the relative strength varies as a function of fre- 
quency, then the spurious RM introduced in the 
measurement could be significant. 

In the right panel of Fig. 5, we have addressed 
the other possibility that for a constant fractional 
strength across the frequency band, (^{v) varies. 
The five curves, solid, dash, dot-dash, dotted and 
dot-dot-dot-dash, correspond to Qiv) varying in the 
ranges of 160°-180°, 165°-180°, 170°-180°, 175°- 
180°, and a constant value of 180°, respectively. 
Here, the RM value changes sign at Fab(z^) = 1, 
and this is exactly what we see in Fig. 2. The 
inferred value of RM((/>) has an anti-symmetry with 
respect to </> ~ 30°, where one mode is stronger 
to the right and the other mode to the left. It 
is worth noting here that in Fig. 5, at the limit of 
FAB{f) — > 00, the RM introduced asymptotically 
reduces to zero. 

6. Discussion 

In Fig. 1 & 2, one of the intriguing things is 

that the value RM(0) does not remain constant be- 
tween various epochs of observation. Perhaps one 
of the reasons could be a combination of interstel- 
lar scintillation and the effect that we have ad- 
dressed in this paper. Although scintillation is 
not expected to have any effect on the amount 
of Faraday rotation introduced in the interstellar 
medium, when the total integration time is short 
(~ a few minutes), the effective centre frequency 
and bandwidth is expected to vary. Combined 
with the effect that we have described, the RM 
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Fig. 5. — Effective RM introduced by non-orthogonal mode- mixing (Y-axis). Left panel shows RM as a 
function of ({v). For a constant C(i^) value, the five curves, solid, dash, dot-dash, dotted and dot-dot-dot- 
dash, correspond to variations of fractional strengths of the two modes (Fab{i^)) from the lower to the upper 
end of the frequency band of 0.3 0.9,0.45 0.9,0.6 0.9,0.75 0.9 and a constant 0.9, respectively. In the right 
panel, the five curves correspond to the variation of C(i^) of 160°-180°, 165°-180°, 170°-180°, 175°-180°, 
constant value of 180°, respectively. See text for details. 



value inferred in principle can be different at differ- 
ent epochs. Moreover, as described in Sec. 2, each 
of the four panels (a to d) in Fig. 1 has been pro- 
duced with 120 second (^200 pulses) long scans. 
Pulsars are known to exhibit stable average-pulse 
profiles only after integrating a few thousands of 
pulses. As we can see from the four panels, av- 
erage profiles, and even polarisation-sweep curves, 
are not identical between all the panels. This may 
also introduce time-dependence in RM values. 

As mentioned in Sec. 5, the exact frequency de- 
pendence of Fab{i^) and ({u) are not known. In 
general, there is no reason to conclude that the 
frequence dependence of (xab(j^)) introduced by 
this effect is the same as that of the RM. What we 
have attempted in Fig. 5 is to estimate the effective 
RM introduced by this effect. A thorough inves- 
tigation of the nature of the emission modes and 
their frequency dependence is beyond the scope 
of this publication, and we plan to consider these 
matters in a subsequent report. 

As radiation from pulsars is highly polarised 
and as distance estimates to pulsars are more reli- 
able than the other polarised sources (e.g. super- 
nova remnants), pulsar RM determinations have 
been used extensively in probing the magnetic field 
structure of the Galaxy. RM is an important tool 



to model the long range as well as the turbulent 
component of the magnetic field in the Galaxy. If 
the effect that we have discussed in this paper is 
common, then it will have serious consequences 
on the existing RM measurements of pulsars, and 
thereby on the models of magnetic field structure 
in the Galaxy. Apart from our detailed analysis 
of PSR B2016-I-28 presented here, in the sample 
of Weisberg et al. (2003), we looked at pulsars 
B0301-hl9, B0525+21, B0626-h24, B1929+10 and 
B2020-I-28. These all showed systematic variation 
of the RM across the pulse profile. For instance, 
Weisberg ct al. measure KM values for these pul- 
sars as -5.7 (±10), -39.2 (±10), 69.5 (±10), -5.9 
(±5), -73.7 (±25) rad m'^, respectively. The val- 
ues given in parenthesis are the approximate mag- 
nitude of variations seen around the mean value 
across the pulse profile. As we can see clearly, 
these variations are significant. Especially, when 
the magnitude of RM is small, these variations can 
cause significant bias. 

The presence of orthogonal modes is common 
among pulsars. Among pulsars for which any 
amount of detailed single-pulse study has been 
done so far, it is clear that a great majority of them 
exhibit orthogonal modes. In particular, it is al- 
most impossible to find a "conal" pulsar that does 
not show evidence of orthogonal mode emission 
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associated with its conal components (Rankin & 
Ramachandran 2003). It is also possible that these 
two modes have significantly different spectral in- 
dices. For instance, there is clear evidence to show 
that the degree of polarisation in the average-pulse 
profiles decreases with increasing frequency. This 
can be understood easily by having one of the 
modes dominating in strength at lower frequen- 
cies, and the two modes having roughly compara- 
ble strengths at very high frequencies. As the ear- 
lier studies (Stinebring et al. 1984; Gil, Snakowski 
& Stinebring 1991; Gil et al. 1992; McKinnon 
2003) have shown, these modes are slightly non- 
orthogonal, and this non-orthogonality of the two 
modes is a wide-spread phenomenon. 

As we have shown in this paper, the only way to 
eliminate this artifact is to determine RM values 
from single pulses. This, of course, is a very chal- 
lenging task, owing to signal-to-noisc; ratio consid- 
erations. For weaker pulsars for which we cannot 
obtain good single pulse data, it is impossible to 
separate the two modes to unambiguously deter- 
mine RM values. Therefore, weaker pulsars are 
bound to suffer from this artifact, and there is no 
obvious way of correcting it. 

On the Galactic scale, whether or not this ef- 
fect will make serious changes to the magnetic field 
model, remains to be seen. A project to deter- 
mine correct RM values for several other pulsars 
is underway, and will be presented in a subsequent 
publication. It is conceivable that this artifact will 
be most prominant among pulsars with small RM 
magnitudes (<50 rad m^^ or so). A thorough 
analysis to check the validity of the already ex- 
isting RM values, and the effect on the Galactic 
magnetic field structure is much needed. 

7. Conclusion 

Our major conclusions from this work can be 
summarised as follows: 

• We find that the rotation measures deter- 
mined for PSR B2016+28 as a function of 
pulse longitude varies significantly. This 
seems to be the case for five other pul- 
sars, namely PSRs B0301-M9, B0525-h21, 
B0626-I-24, B1929-M0 and B2020-I-28. 

• This effect is an artifact introduced by the 
frequency dependence of relative strengths 



of the two modes as well as the amount of 
non-orthogonality, and it is not intrinsic to 
the pulsar magnetosphere. We show that 
if we estimate rotation measure for the two 
modes separately, this effect can be removed. 
The technique to remove this artifact in- 
variably involves analysis with single pulses. 
Hence, cannot be carried out for fainter ob- 
jects. There is no other obvious way of com- 
pensating for this artifact for these fainter 
objects. 

• As the amount of RM spuriously introduced 
can be as high as a few tens of units, several 
measurements of Rotation Measures of pul- 
sars in the literature may be in error, and 
needs revision. 

We thank Jon Arons, Avinash Deshpande, Aris 
KarastergioU; Simon .Johnston, Michael Kramer 
& Dipanjan Mitra for useful comments on the 
manuscript. Special thanks are due to Dipanjan 
Mitra, Simon Johnston and Michael Kramer for 
sharing their results prior to publication. We also 
thank our anonymous referee for his/her valuable 
comments. 
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DEFINITIONS 







1") p fi n 1 i", 1 o n 

1 ^ %^±±±±± l,;X\^±± 




Y 






PA 






RM 


l-?nf fif inn A/ToaQiirp 






r?n1"fii'if~m \l\ c;n fa fi i" fi cri'\7'A'n 1 f~m ccTfii Hp /ti tti as qi i vpH n ?i'\7'PTfi cp tii i 1 cp T^vr^ri Ipq 

1\,(J UdtlUil lVJ.CCLOU.iC dt d ^IVCli iUil^i t LltJ-C iliCdO Lli CU. Witli dVCid^C JJLlloC ^HJillCiD. 




nMA{<l>) 


Rotation Measure of polarisation mode-A 




UMbW 


Rotation Measure of polarisation mode-B 




SNR 


signal-to-noise ratio 


XA{(t>, 




Weighted mean values of linear polarisation position angles of 






Mode-A and Mode-B 




(xab(j^)) 


linear polarisation position angle at a given frequency after 






incoherently superposing the two (non-)orthogonal modes. 


XA{t, 




linear polarisation intensities of the two modes 






as a function of time and frequency. 




X^it, v) 


system noise power as a function of time and frequency. 




(RM) 


Net RM value computed from average pulse profile after averaging 






over all pulse longitudes. 
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